arXiv:1508.04452v2 [quant-ph] 3 Feb 2016 
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For N spin-polarized fermions with harmonic pair interactions in a 1-dimensional trap an odd- 
even effect is found. The spectrum of the 1-particle reduced density matrix of the system’s ground 
state differs qualitatively for N odd and N even. This effect does only occur for strong attractive 
and repulsive interactions. Since it does not exists for bosons, it must originate from the repulsive 
nature implied by the fermionic exchange statistics. In contrast to the spectrum, the 1-particle 
density and correlation function for strong attractive interactions do not show any sensitivity on the 
number parity. This also suggests that reduced-density-matrix-functional theory has a more subtle 
A^-dependency than density functional theory. 
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Introduction .— Physical behavior often depends 
qualitatively on binary parameters as, e.g., odd/even or 
integer/half-integer. Such parity effects play an impor¬ 
tant role in physics. A well-known example is Kramers’ 
number-parity effect, i.e. the twofold degeneracy of the 
eigenstates of a quantum system with an odd number 
of electrons, provided time reversal symmetry holds [I|. 
Haldane Q has shown the existence of a spin-parity ef¬ 
fect. The spectrum of the quantum Heisenberg antifer- 
romagnet in one dimension has an energy gap for all in¬ 
teger spins whereas it is gapless in case of half integer 
spins. Recently, an interesting number-parity effect has 
been observed experimentally for a few ultracold fermions 
in a quasi-one-dimensional trap. Tuning the potential 
such that the pair interactions become attractive, Cooper 
pairs are formed. Their tunnelin g i s different for an odd 
and an even number of fermions |3|. Based on Kramers’ 
theorem another number-parity effect was proven to ex¬ 
ist for fermionic 1-particle reduced density matrices (1- 
RDM) Q (see also Ref. 0)- There it was shown that the 
eigenvalues of a 1-RDM arising from an eigenstate of a 
time reversal symmetric hamiltonian are twofold degen¬ 
erate for an even number of fermions. 


In the present work we will show that the so-called nat¬ 
ural occupation numbers, i.e., the spectrum {Xk} of the 
ground state 1-RDM for strongly attractive and spinpo- 
larized fermions confined in one dimension also exhibits 
an odd-even-effect. Since the spin-polarizing magnetic 
field breaks time reversal symmetry, this effect is com¬ 
pletely different from that found in Ref. Q . Furthermore, 
it does not occur for bosons. Consequently, it must result 
from the fermionic exchange symmetry. 


Besides the relevance of parity effects on their own, ex¬ 
ploring the structure of reduced density matrices has also 
gained a lot of relevance during recent years. This is es¬ 
sentially due to progress MU in the quantum marginal 
problem (QMP) which studies the relation of reduced 
density matrices arising from a common multipartite 
quantum state. For basic overviews of the QMP the 
reader may consult Refs. 12h14|. The most prominent 


QMP is the 2-particle A’-representability problem BQ, 
the description of 2-RDM arising from A-fermion quan¬ 
tum states. Its solution would allow to efficiently calcu¬ 
late ground states of fermionic quantum systems with 2- 
body interactions. However, since this problem and most 
of the other QMP are Quantum-Merlin-Arthur-hard [l^ , 
already partial insights on the set of compatible density 
matrices are highly appreciated and alternative methods 
for the ground state calculation are gaining importance 
as well. One such promising method is reduced-density- 
matrix-functional theory (see e.g. |l7l IlSl l. This natural 
extension of density functional theory [l9| seeks a distin¬ 
guished functional J- on the 1-RDM whose minimization 
leads to the exact ground state energy and the corre¬ 
sponding ground state 1-RDM. Any structural insights 
on ground state 1-RDM contributes to this task of find¬ 
ing or approximating A by exposing further necessary 
constraints on legitimate functionals. 

Model and 1-RDM .— We consider A identical par¬ 
ticles with mass m in a one-dimensional harmonic trap 
interacting via a harmonic two-body potential. If Xi is 
the position of the Tth particle the hamiltonian reads: 


N 




i=l 


1 2 2 
2m dx^ 2 




E 


Xi-X, 


l<i<j<N 


( 1 ) 

where u is the eigenfrequency of the trap and D the inter¬ 
action strength, which can be positive or negative. Sta¬ 
bility requires D > Di^w = —moS^/N. 

Hamiltonian © arises as an effective model, e.g., for 
the description of quantum dots, where the Coulomb 
interaction between the electrons is screened (see 
e.g. Ref. 0 ). Furthermore, it was used to under¬ 
stand the emergence of shell structures in atoms (see 
e.g. Ref. [2l[) and nuclei (see e.g. Ref. (22l|b 

The great advantage of model o is the exact knowl¬ 
edge of all its eigenstates 231427j . For arbitrary numbers 
of bosons and any spatial dimension the ground state 1- 
RDM can easily be calculated HQ. Based on such 
analytic results Bose-Einstein condensation in harmonic 
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traps was explored (see, e.g., Ref. (28|). In contrast to 
bosons, the analytical calculation of the corresponding 
fermionic 1-RDM is much more involved. We will show 
that the properties of the fermionic 1-RDM are much 
richer, compared to the bosonic case, leading to new in¬ 
sights. 

For spin-polarized (or spinless) fermions in one dimen¬ 
sion, the spatial part of the ground state of hamiltonian 
o is the totally antisymmetric wave function |23M25j| 
(see also Refs. EflUl) 


l<i<j<N 

■ exp[— A[xi x‘^) B[xi X]\f)‘^] (2) 


with a normalization factor and 


1 / 1 




2/V ^ 2iVV/2 /2 


(3) 


/_ = ^Jh/mu and /+ = -\- DNlm)^/"^ are the 

length scales for the center of mass and the relative mo¬ 
tion, respectively. Note that ^'g'^^(a;i,..., Xn) resembles 
Laughlin’s wave function [2^ for the fractional quantum 
Hall effect. We will come back to t his p oint below. 

The 1-RDM of ((2]) follows as 24. 25| 


Pi^\x; y) = J\f^-^'^F{x, y) exp[—a(a;^ + y"^) + hxy] (4) 


with the polynomial F{x, y) in x and y of degree 2(iV — 1) 
(see Appendix!^ and 


a = A-B-hl2 , h = {N-l)B'^/[A-{N-l)B]. {F) 


The factor (not to be confused with the nor- 

malization constant of 4/g '^), follows from the normal- 
oo f n 

ization f dxpY (x; x) = N. The bosonic and fermionic 

— OO _ 

1-RDM are related by 


p[-^^ {x]y) = M^f'^F{x, y) p^^^ {x;y), (6) 

where JFB) = J\fB) and the corresponding 

normalization constant appearing in p^^\x;y). Conse- 
quently, the fermionic nature of p\ {x;y) is only con¬ 
tained in the polynomial prefactor F{x, y). R arises from 
the polynomial prefactor of the exponential function in 
Eq. ((2j, the Vandermonde determinant, which is a result 
of the fermionic exchange symmetry. 

The spectrum {Afc} follows from solving the eigenvalue 
equation 


j dypY\x;y)Xk\y) = Xk\k\x), (7) 


A; = 1, 2, 3,.... In quantum chemistry the Afc are called 
natural occupation numbers. These eigenvalues will be 
ordered decreasingly, i.e., Afc > Afc+i for all A: > 1. Due to 
the duality {Afe(/_,/+)} = {Afc(/-i-,/-)}, first observed in 
Ref. [13 and proven in Ref. [3l|, we restrict to /_|_//_ < 1. 
Notice, /+//- < 1 (/+//- > 1) means attractive (repul¬ 
sive) pair interactions. Since we are interested only in 
fermions we also suppress the superscript (/). 

Strong coupling limit and results. — For attractive 
interaction the strong coupling limit DN/muj'^ —> oc 
corresponds to the limit t = 1+/1- 0, performed 

at N fixed. For repulsive interactions it follows for 
DN/moo'^ DiowN/muj^, corresponding to t ^ oo. 

We will prove that completely unexpected features of 
the spectrum of the fermionic 1-RDM occur in this limit, 
which can be discussed analytically. Due to the duality 
property |3l| we can restrict to attractive interactions, 
i.e. to A < 1. Technical details can be found in the 
appendix. 

In the following all coordinates Xi and all lengths will 
be measured in units of /_. For A ^ 0, there are three 
basic observations. First, the leading order of pi{x;y) 

is proportional to exp — ^^^^((x — ?/)/A)^j (see Ap¬ 
pendix [Cj. Therefore, the weight of the deviation of y 
from X decreases extremely fast for decreasing A. Sec¬ 
ond, Xk{x) varies on an x-scale proportional to \/A and 
third, the polynomial prefactor F{x,y) converges to a 
polynomial F{z) where z = {x — y)/t (ses Appendix [Bj. 
Let us introduce the rescaled variable x = xj^Jt and 
XkiVix) = exp(iA’A5^)Cfc(i). Then, by using the mo¬ 
mentum representation the eigenvalue equation (j7]) for 
A ^ 1 reduces to a Schrodinger equation (with position 
and momentum exchanged) 


- I Cfe(P) = i-^k) d(p) . (8) 


Here Cfc(p) is the Fourier transform of C,k{x), 

fhN = h^/{-2tNVN{0)) (9) 


the mass of the ‘particle’ and 

OO 

Vn{p) = —tj\f J dzF{z) cos{'/ipz) exp ^- 

— OO 

( 10 ) 

an effective potential, where the scaled momentum p = 
y/tp/h has been introduced, p is the conjugate momen¬ 
tum of X. Vn{p) is a product of a polynomial in {\/ip)‘^ 
of degree {N — 1) (originating from the fermionic nature 
of the 1-RDM) and a Gaussian ey^p[—j^^{\/ip)‘^] (see 
Appendix 

Figure [1] depicts Vn{p) for V = 3 and N = 4. One ob¬ 
serves that its number of minima equals N. Since V/v(p) 
for arbitrary N is symmetric, one of the minima must be 
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FIG. 1. (Color online) Scaled effective potential as a function 
of Vtp for = 3 (left panel) and = 4 (right panel). 

at p = 0, for N odd. In this case it is the global mini¬ 
mum, which we checked systematically up to A?” = 19 and 
for some exemplary up to = 101. For N even, the 
global minimum is twofold degenerate, which we checked 
for N = 2,4, 6,..., 20 and for some exemplary N up to 
N = 100. There is little doubt that these properties hold 
for all N. This qualitatively different behavior of Vn{P) 
implies that the spectrum {A/.} of the ground state 1- 
RDM will qualitatively differ for an odd and even num¬ 
ber of particles. Without solving the eigenvalue equation 
(HI) one can already predict that the spectrum a of the 
1-RDM consists of two parts for k < and 

cr(^oui) ^ ^ k^{t) where k^{t) 1/t is the A;-value for 

which (—Afc) equals the height of the highest maximum 
of Vn{p)- CTen^^ consists of “isolated” eigenvalues, only, 
and does not depend qualitatively on the number parity. 
In contrast, the upper part, o'^\ differs qualitatively for 
N odd and even. For N odd it consists of subsets with 
“isolated” eigenvalues and subsets of pairs of quasidegen¬ 
erate eigenvalues, whereas for N even only subsets with 
quasidegenerate eigenvalues occur. 

This number parity effect can be illustrated by calcu¬ 
lating analytically the largest eigenvalues of the 1-RDM 
corresponding to the low-lying eigenvalues (—Afc) of the 
‘particle’ in the effective potential. This will be done by 
use of the harmonic approximation Vn{P) — + 

— pminY for the global minimum. Then, 
Eq. (HJ reduces to the eigenvalue equation of a harmonic 

oscillator with mass m, frequency O = /fh and 

eigenvalues Afc = — (Afc + = hQ{k — Accord¬ 

ingly, we obtain 

Afc ~ t ajl - t( 11 ) 
for A: = 1, 2, 3,... . The coefficients follow from 

a = /3 = , ( 12 ) 

which are of 0{t^). 

For N odd, the global minimum of Vn (p) is nondegen¬ 
erate at p = 0 . Hence the smallest eigenvalues (—Afc) and 
therefore the largest eigenvalues Afc of Eq. ((7|) are “iso¬ 
lated” . In case of N even, the global minimum is twofold 


degenerate such that the tunneling between both min¬ 
ima leads to a splitting (AA^ + AA^). Therefore, the 
eigenvalues appear in quasidegenerate pairs {A 2 fc-i, A 2 fc} 
given by 

Asfc-i ~Aa[l-A/3(/c-i)] +AA+ (13) 

for A; = 1, 2, 3,... . X 2 k follows from Eq. (fT3l) by replacing 
AA+ by (-AAfc). 

The results OD and o show that Afc is of 0{t). Its 
^-dependency is proportional to A^. Since these eigenval¬ 
ues are the occupancies of the 1 -particle states Xk{x) they 
are nonnegative. Hence, the validity of Eqs. HU), HI is 
limited to k < k^ ^ 1/t. They are also based on the har¬ 
monic approximation, which restricts their validity even 
more(see below). Furthermore, the tunneling splitting 
(AA+ + AXp) can be estimated. Equation OT shows 
that V]\f{p) = 0{t) and that it varies on a scale l/s/i- 
This implies m = a potential barrier AVn = 

0{t) and a tunneling distance Ap = 0{l/\/t). Using that 
the splitting is proportional to exp {—y/rh AVn Ap) / h) h 
follows that (AA^ + AA^) ^ exp (—0(1/A)), i.e. for N 
even the pairs (A 2 fc-i,A 2 fc) become quasidegenerate in 
the regime of strong coupling. 

In order to test these predictions we have deter¬ 
mined the eigenvalues by solving numerically the orig¬ 
inal eigenvalue equation o for AI = 3 up to A = 8 . 
Note that N = 2 is a special case, since there, all 
eigenvalues are automatically twofold degenerate as a 
hamiltonian-independent consequence of the fermionic 
exchange statistics (see Theorem 4.1. in Ref. (H). Figure 
[Hlpresents the larger eigenvalues for the coupling t = 10“^ 
for A = 3 and A = 4. In case of A = 3 the subsequent 
eigenvalues have almost the same distance which is of 
O(A^) whereas for A = 4 they occur in quasidegener¬ 
ate pairs with A 2 fc-i - A 2 fc+i = O(A^) and A 2 fc-i - X 2 k ^ 
exp (—0(1/A)), as predicted by our analysis above. Com¬ 
paring for A = 3 the numerical result with the analyti¬ 
cal one, Eq. CD, shows very good agreement for Ai and 
the slope dXk/dk, for small k. For A = 4, Ai is also 
well reproduced. However, the analytical and numerical 
values for the slope dXk/dk deviate stronger from each 
other. This is a consequence of the fact that in contrast 
to A = 3 the harmonic approximation is rather poor due 
to (p — prnm)^-anharmonicities of U 4 around pmm- Note, 
our major achievement is not result m and 03 for the 
largest eigenvalues, but 

(i) the qualitatively different spectrum of the 1-RDM 
for A odd and even in the strong coupling limit 

(ii) the generation and annihilation of quasidegenerate 
pairs of eigenvalues for those indices k for which 
(—Afc) becomes equal to the height of corresponding 
minima and maxima, respectively, of lAv(p) 

Figured] demonstrates the creation of pairs of quaside¬ 
generate eigenvalues and its inset illustrates their annihi- 








4 




FIG. 4. (Color online) x- and y-dependency of the 2-particle 
density n{x,y) (not normalized and we set I- = 1) for t = 0.1, 
N = 3 (left panel) and = 4 (right panel). 



k 

FIG. 2. (Color online) Spectrum {Afe} for t = 10^^. Upper 
panel: N = 3 and I < k < 40, lower panel: N = A 

and 1 < A: < 62. The solid lines present the correspond¬ 
ing analytical result for {\k} from Eq. dill) and from the 
first term of Eq. m, respectively. The insets show the h- 
nal crossover from the regime of quasidegenerate to that of 
non-quasidegenerate eigenvalues. 
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FIG. 3. (Color online) x-dependency of the 1-particle density 
for = 3 (left panel) and N = A (right panel), where we 
set I- = 1. Blue dashed, black solid and red dotted line 
correspond to t = 0.1 (attractive), t = 1 (free fermions), 
and t = 10 (repulsive), respectively (for f = 10 it is plotted 
lOn(lOx)). 


lation, e.g. at the highest maximum of Vn{p), where the 
crossover to “isolated” eigenvalues of occurs. 

Since under an increase of N, Vn{p) develops more and 
more extrema, there will be more and more regimes with 
groups of a different number of quasidegenerate pairs. In 
order to resolve these regimes for large N, t must become 
small enough. Since the t and N dependence occurs as 
tN (see Appendix IXI) it must he t t^{N) ^ 1/N. For 
macroscopically large N of 0(10^^) this requires values 


for t which are not realizable in experiments. Yet, for 
N of 0(10^), for which already macroscopic properties 
are present, this should be feasible. Accordingly, few- 
fermion systems are particularly suitable to observe these 
qualitative features of the spectrum of the 1-RDM. 

\k ior k ^ ^ 1/t can be determined analytically. 

With the approach discussed in Ref. [2^ we obtain (see 
Appendix [Cj 


Afc -- {kt) 


N-l 


r 2Y 1, 

exp- , t [k — -, 


( 14 ) 


independent of the parity of N. 

We have also calculated analytically the 1-particle 
and 2-particle density n{x) = pi{x]x) and n{x,y) = 
p 2 {x,y]X,y), respectively. The latter follows from the 2- 
RDM p 2 {x, y; x', y'). n{x) is shown in Figure[3]for N = 3 
and N = A. For noninteracting fermions (black solid 
line), the ‘layering’ of the particles within the harmonic 
trap can be seen. This “shell structure” also exists for 
repulsive (red dotted line) and weak attractive coupling, 
but disappears completely for strong attractive interac¬ 
tions (blue dashed line), becoming qualitatively indepen¬ 
dent of N, which we checked up to Y = 20. The duality 
discussed in Ref. [3l| implies that the 1-particle density 
in momentum space for strong repulsive interaction be¬ 
haves similarly as n{x) for strong attractive interactions, 
i.e. it becomes structureless, as well. Quite similar be¬ 
havior has been found for n{x,y). As demonstrated by 
FigurelUthe ‘layering’(not shown) for n{x,y) disappears 
again for strong attractive coupling, and does not ex¬ 
hibits any qualitative sensitivity on N (cf. left and right 
panel of Figure |4|). All these properties also hold for the 
correlation function C{x,y) = n{x,y)/n{x)n{y). 

Summary and Conclusions .— We have shown that 
the 1-particle description in form of the 1-RDM exhibits 
an odd-even effect. The spectrum of the 1-RDM related 
to the fermionic ground state of our one-dimensional har¬ 
monic system differs qualitatively for an even and an odd 
number of particles. This effect does only occur for strong 
attractive and repulsive (due to the duality property (SlJ) 
interactions. The number-parity effect does not exist for 
bosons. Therefore it must originate from the repulsive 
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nature implied by the fermionic exchange statistics (an¬ 
tisymmetry of the wave function). One may wonder how 
far the interplay between strong pair interactions (par¬ 
ticularly for attractive ones) and the exchange symmetry 
leads to new phenomena for fermions, beyond the present 
parity effect. Also the investigation of the existence of the 
odd-even effect in more than one spatial dimension will 
be of interest. 

It would also be interesting to develop tools which 
make it possible to investigate these predictions by ex¬ 
periments. In that case the parameter t has to be tuned 
(see below) such that the splitting of the quasidegenerate 
eigenvalues is still large enough in order to be resolved. 
Since our model involves harmonic pair interactions one 
might be tempted to deny the relevance of our findings 
for realistic systems. There are two reasons why this 
might be not true. First, for arbitrary pair potential the 
particles in a trap form a one-dimensional lattice. Ex¬ 
panding the potential up to quadratic terms in the dis¬ 
placements with respect to the classical groundstate will 
result in a harmonic model similar to that studied by 
us. This has recently been done for a one-dimensional 
A^-particle system with long-range inverse power-law po¬ 
tential in order to calculate the von Neumann entangle¬ 
ment entropy 3^. Second, and even more important, the 
specific form of the pair interactions may not be as im¬ 
portant as one might believe. As pointed out above the 
odd-even effect for the spectrum of the fermionic 1-RDM 
is a result of the polynomial prefactor in Eq. 0 which 
makes the wave function totally antisymmetric. This re¬ 
quirement of antisymmetry has also been the guide lead¬ 
ing to Laughlin’s wave function for describing the frac¬ 
tional quantum Hall effect. This wave function again 

with 


[<i<j<N 


{Zi-Zj)Pj, 


has a preexponential factor 
p an odd natural number and Zi the complex variable 
specifying the position of the i-th particle in the plane. 
Laughlin’s ansatz is a surprisingly good approximation 
of the two-dimensional electron system’s ground state, 
not only for a Coulombic pai r potential but for har¬ 
monic interactions, as well |29|. This suggests that the 
specific form of the pair interactions is less important 
which is supported by the exact ground state solution 
of the one-dimensional Calogero-Sutherland model 33 


37l| . Besides a harmonic trap potential this model con¬ 
tains pair interactions g{xi — Xj)~^. Its ground state 




involves a preexponential factor 

with p{g) = i(l ± + 4g). Choosing the coupling con¬ 

stant such that p{g) is an odd natural number one obtains 
the ground state for N fermions. 

Therefore our results may also hold for ultracold 
fermionic atoms in an optical trap interacting by a con¬ 
tact potential . The interaction can be tuned from 

attractive to repulsive. In particular it can be made ar¬ 
bitrary strong by approaching the Eeshbach resonance. 
This would allow to experimentally approach the strong 


coupling limit discussed in the present paper . Of course, 
another possibility to realize that limit is the decrease 
of the trap frequency. Eor instance for Li the trap fre¬ 
quency in the experimental setup in Ref. |^lis about 
thirty times smaller than in the setup of Ref. ^ . 

The odd-even effect may also initiate and guide a 
new direction in density and reduced-density-matrix- 
functional theory. Although three decades ago, it had 
been argued that the A’-dependency of the ground state 
energy implies an A-dependency of the density function¬ 
als Aiv |43|, all of the prominent functionals used today 
do not exhibit an explicit dependency on A (see e.g. (3). 
The parity effect found by us is the first demonstration 
of a subtle A-dependency of the spectrum of the ground 
state’s 1-RDM, and therefore also of TA. In addition, 
the fact that the 1-particle density n{x) in the strong cou¬ 
pling regime does not show any sensitivity on the number 
parity suggests that the A dependence within reduced- 
density-matrix-functional theory [13, is much more 
subtle than in ordinary density-functional theory [l3 |. 
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Appendix A: Calculation of F{x,y) and Vn{p) for t ^ 0 


In this appendix we calculate the leading order in t of the polynomial F{x^ y) and of the effective potential Vn{p)- 
F{x, y) is given by 


N-l 


F{x,y)= I due ^^Hk{pu + q{x,y))Hk{pu + q{y,x)), 


k=0 


where Hk{x) are the Hermite polynomials and 

p = [B/{A -{N - , q{x, y) = - ]^P^{x + ?/)j . 

From Eq. ()A2p we obtain with Eq. o and t = /_(_//_ 


1 


p = 1 + 0{t ), q{x, ?/) = - 1 + 0{t ) —--h 0{t)x . 


2A^-y 


t 


The following identity 


Hk{x + y) = (^^Hm{x){2y) 


m=0 


\k—m 


(Al) 


(A2) 


(A3) 


(A4) 


is useful. It can be proved by Taylor-expanding its left hand side with respect to y and using for the n-th derivative 


H^^\x) = 2'^jB^^Hk-n{.x) which is easily proved by induction taking the initial condition = 2kHk-i{x) 


fc! 


(fc—n)! 




45l | into account. Substituting the leading order term of p and q{x^ y) into Eq. (|Aip and making use of the identity 
one gets 


N—1 ^ k k 

F{x,y)^J2F^Yl 


2^k\ 

k=0 m=0 m'=0 


^ \ f ^ \ f ^ ~ fy ~ X'' k-m' 


t 

du C B j-fiiF) Byyi! {u) 


N-l ^ k 
k=0 m=0 


^Y2rn^l( 
mJ V 


N-l 


= Vtt 


- 1 ) 


iV-l 


2i^m\ I ^ 


m=0 

N-l 

= ^/?E 

m=0 

^ m , 

CX) 


k\l/X — y 


m 


_ 9 ; \ 2m 


t 


k=m 

(_l)m ^ 


2™-m! Vm + 17 V t 


where 5 = {x — y)/t. We have used / du 6““ Bm{u)Bm'{u) = y/F2'^m\5mm' 451 and ^ ^ 01 • Note 


(A5) 


that F{z) is an even function in 5. It is straightforward to calculate {z^'^F{z)) = J dzz^'^F{z) exp One 

— OO 

obtains 

N-l 

(V(r) -I- rr) 1 — MM / /V \ / V/V \ m+n 

(A6) 


m=0 


Making use of Eq. ()A5p one can also calculate the effective potential (Eq. (fTU|) f which can be represented as Vm{p) = 
(F(5) cos(-\/tp5) exp( —One obtains with A/” ~ \/]V/7r(cf. Eq. 


N-l 


Vn{p) = -^[ 


tN 


(A7) 


n=0 





















with coefficients 


'<^N,n 


2iV (2(m-n)-1)!! / N w2m\/ 2N 

- 1 ^ ^ 2"^m! Vm + 1A 2n AiV - 1/ 

v m=n 


(A8) 


Appendix B: Mapping of the 1-RDM for t ^ 0 


This appendix presents details of the reduction of the eigenvalue equation ([7]) for the 1-RDM in the strong coupling 
limit t ^ 0 to an eigenvalue equation for a particle in an effective potential Vn{p)- Using Eqs. ((3]) and ([5]), one obtains 
for the exponent on the right hand side of Eq. o 


a{x^ + ?/^) - bxy = (^1 + O(t^) 


x-y 

t 


+ N{1 + 0{t'^))xy . 


Erom Eq 


we have 


N-l 


F{x,y) ^ 


m=0 


2"^m! Vm + 1 


X - ?/\2m ^ - y 


(Bl) 


(B2) 


Taking into account that Eqs. (|B1D and (|B2|) imply a(x^ + — bxy]x=y = Nx"^ and F{x,x) = \FtN^ respectively, 

the normalization condition for pi{x; y) gives the normalization constant in leading order in t 


Af ^ 


^/N 


(B3) 


The results CD and (|B2I) show that pi{x; y) in leading order only depends on {x — y), i.e. is translationally invariant. 
This can easily be understood since for /_|_ fixed, t ^ 0 implies that the length scale /_ associated with the harmonic 
trap converges to infinity, i.e. its eigenfrequency (u goes to zero. Consequently, the external potential becomes fiat for 
t ^ 0, and the hamiltonian m becomes invariant under translations. To proceed, we introduce scaled variables. 


Then we get from Eqs. 0, 


= x/Vi , y = y/Vi , z = {x-y)/t = {x-y)/Vi. 


N-l 
4N ' 


pi(x;y) = AAexp (—exp (iVt2x5)E(5) exp — 


(B4) 


(B5) 


Going beyond the leading order term (|B2p of F{x,y) there also appear additional terms {{x — y)/t)'^ = z'^ for 
m = 0,1,...,A^ — 1 with coefficients of order and terms x{x — y/tz) with a coefficient of order Therefore the 
right hand side of Eq. ()B5I) with F{x,y) given by Eq. (|B2I) represents the leading order terms of the expansion of 
/9i(x, y) for t ^ 0. 

Next, we introduce scaled natural orbitals 


By making use of Eqs. 
(iV - l)/(4iV)) =a 


Xk{x) = XkiVix) . (B6) 

and choosing 5 as an integration variable the eigenvalue equation 0 becomes with 


t Af / dzexp{—Ntx‘^)exp{NFxz)F{z)e Xk{x — Viz) = \kXk{x). 


(B7) 


The second exponential in Eq. 


This leads to 


with the exponent involving xz can be eliminated by using 
»(i) = exp (lNtx^)(k{x) . 


(B8) 


t N'expi—Ntx^) / dzF{z)e Cik{x —'/tz) = \kCk{x) ^ 


(B9) 












' 1 AT-|-z2^ 

Next we use 


where exp {^Ntz'^) (which occurs in the integrand) can be neglected with respect to exp (—since a = 0{t^). 


Substituting Eq. ()B10|) into Eq. 


operator v = 

^ ^ % OX 


Cfc(5 - V^5) = exp ^Cfc(^) • (BIO) 

and taking into account that F{z) is even in 5 we get with the momentum 


t Af exp{—Ntx‘^)^ f dzF{z)e cos {Vizp)^Ck{x) = \kCk{x) 


(Bll) 


Using momentum representation, i.e. Eourier transforming Cfc(5) 


Cfc(i) = ^ 


(B12) 


turns out to be the appropriate choice. Then Eq. (IBllI) becomes 

CXD 

/ d5.P(5)cos(Vt5p)e““"'''JCfc(p) = (-A/c)Cfe(p) • 


exp I 


(B13) 


Note, the term in the square bracket is just the effective potential Vn{p) from Eq. (fTOl) . The final step is the expansion 
of exp {Nt-§^). This yields exp {Nt-^)VN{p) = [1 + + 0(t^)]Uiv(p) ~ Vn{p) + Vn{p = 0)iVt^. Substituting 

this into Eq. ()B13I) leads immediately to the eigenvalue equation ([8]) with the mass m given by Eq. ([9]). 

Einally, we justify the scaling ()B6P and the validity of the implicit assumption that Xk varies on a length scale 

X = 0(1). Erom Eq. fTO|) one gets U/v(0) = 0(t), = O(t^) which leads to m = 0(t“^) and O = = 

0{t^). The length scale corresponding to the harmonic oscillator equation obtained from Eq. (j8]) by expanding Vn{p) 
up to second order in {p — Pmin) is given by L = w which is 0(1). Accordingly, (k, and therefore 
scale of 0(1). 


Xk, varies on a 


Appendix C: Calculation of Xk for k ^ k*{t) 


Here we follow the strategy used in Ref. [2^. We present the crucial steps, only. Eor details the 
Appendix C of Ref. [2^. There it was shown that pi{x; y) = {x\pi\y) with the 1-particle reduced 


reader may consult 
density operator 


N-l 2v 

pi = 

1>=Q ij.=0 




(Cl) 


and the effective hamiltonian HeS = at ^ 0^0 + a and are boson annihilation and creation operators and 

the coefficients of the polynomial F{x, y). The explicit form for Oat is given in Ref. 2^ and is not essential here. The 
position operator x (in units of /_) can be represented as follows 


X = —^ (a + a^) = + o^) 


V2 

with the bosonic length scale and = {N — 1 )~^F _ 

Let |m) be an eigenstate of a^a with eigenvalue m and m 1. Then it follows for p <C m 


(C2) 


H- 

(a + a'*')^|m) = |m + p — 2 k) . 


K = 0 


(C3) 
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Following Appendix C of Ref. 


25l | and correcting a few typos one obtains 



r=-{N-l) 


where m ^ t ^ has been assumed, for fixed N. The coefficients {hm,m+r} can be expressed by 
depend explicitly on m. Therefore, the eigenvalue equation 


PilXfc) = Afcixfc) 


for \xk) = Xk^'^l'iTT') becomes a finite difference equation for the coefficient Xk 


(m) 


m=0 


r=-[N-l) 


In Ref. [2^ the eigenvalues \k in Eq. ()C5|) were determined. Using the leading order result — 

(which follows from Ref. [IBj) one obtains 




(C4) 

and do not 

(C5) 
2 iV(Z^))2 t 

(C6) 


which is up to a constant identical to Eq. aa. 



